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Abstract 

We report on the closing of gaps in the ground state of the critical Heisenberg- 

Ising chain at momentum vr. For half-filling, the gap closes at special values 

of the anisotropy A = cos(7r/Q), Q integer. We explain this behavior with 

the help of the Bethe Ansatz and show that the gap scales as a power of 

the system size with variable exponent depending on A. We use a finite-size 

analysis to calculate this exponent in the critical region, supplemented by 

perturbation theory at A ~ 0. For rational 1/r fillings, the gap is shown to 

be closed for all values of A and the corresponding perturbation expansion in 

A shows a remarkable cancellation of various diagrams. 
72.15Nj, 05.30.-d, 75.30Ds 
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I. INTRODUCTION 



In recent years, there has been much interest in the finite-size spectra of one-dimensional 
(ID) quantum many-body Hamiltonians. This effort was largely motivated by the fasci- 
nating ideas of conformal invariance, which allow for a computation of scaling indices and 
classification of universality classes directly from the 1/ L behavior of a given critical model 

i- 

For ID systems, solvable by the Bethe-Ansatz (BA) method, the finite size spectra are 
accessible analytically by various methods and these models have thus been used as testing 
grounds for the general ideas of conformal invariance. In particular, the Heisenberg-Ising 
(H-I) model of an anisotropic spin-1/2 chain has been studied in great detail It has 
been shown that its low energy spectrum in the large N limit can be classified as that of a 
c = 1 conformal field theory with marginally irrelevant operators giving rise to logarithmic 
corrections [0]. 

In the present work, we shall reexamine the behavior of the spectrum and the correspond- 
ing states of the H-I chain threaded by a magnetic fiux $. As shown in 0, the adiabatic 
ground state for half-filling is periodic in $ with a period equal to 2 ■ 27r — 27r is the fiux 
quantum in appropriate units — for values of the interaction strength A = cos{tx/Q) with 
Q not an integer. This is due to the presence of a finite gap at $ = 27r. We show that 
this gap scales as a power of the system size, with variable exponent depending on A. For 
A ~ 0, we calculate the exponent explicitly by perturbation theory. For Q equal to an 
integer, however, the gap closes exactly and the periodicity of the ground state becomes 
macroscopic at half-filling, i. e. proportional to the system size. We show how this behavior 
can be understood in terms of the BA solution as a successive crossing of string solutions. 

Furthermore, for 1/r filling with r an integer, the periodicity is again not the expected 
r fiux quanta and we observe exact degeneracies between ground state and first excited 
state at various values of $. This is true for all values of the interaction strength and not 
only at isolated points. We show that the numerical results are confirmed by second order 



2 



perturbation theory, due to a remarkable cancellation of various diagrams. 

Lastly, we demonstrate that all degeneracies can be removed by simply adding a non 
integrable next-nearest neighbor term to the Hamiltonian. 

II. CLOSING OF ENERGY GAPS AT HALF-FILLING 

As is well known, the H-I model has an interpretation as a lattice gas of either fermions 
or bosons, where spin-down represents a particle and spin-up represents an empty site. The 
magnetic flux $ modifies only the hopping term, s.t. the Hamiltonian is 

N N 

H = -- ^{e^^/^cjcj+i + h.c.} + cos/i ^ ^j^j+i + eo, (1) 

^ 3=1 3=1 

where Cq = A(2M — N/2), N is the number of sites, M the number of particles and we 
reparametrized the anisotropy as A = — cos(yu). We note that the total momentum induced 
into the system by the flux is given as P = M^/N. We further remark, that the repulsive 
region < /i < 7r/2 and the attractive region 7r/2</i<7rof the H-I chain are related by 
an inversion of the spectrum: The behavior of the low-lying energy states in the attractive 
region will be mirrored to the behavior of the high energy states in the repulsive region and 
vice versa. 

Let Eq{<^) denote the ground state energy of the system with flux $. As is easy to see, 
i?o($) has a periodicity of 27r in $, i.e. Eq{(^ + 27r) = i?o(-$) = Eo{(^). However, the 
periodicity of a given state can be an integer multiple of the period of the energy. Following 
0], we define a topological winding number n to be the number of times the flux $ increases 
by 27r before the state returns to its initial value. 

We write the shift of the ground state energy as a function of flux as A£^o(^) = Eq{(^) — 
£;o(0) = D^'^/2N + 0(<l>^), where D has been called the stiffness constant §]. For < /i < 
7r/2 and a half-filled band (M = N/2), the stiffness has been calculated exactly as = 
Vs/2{iT — fi) for flux values |$| < 2(7r — yu); the spin- wave velocity is given by f ^ = vr sin fi/jj, 

0- 
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At half-filling, we now adiabatically boost the zero momentum ground state until we 
have $ = 27r and P = n. For < /i < 7r/2, we then observe a finite gap AE between 
the boosted ground state and the first excited state in this momentum vr sector and so the 
winding number of the zero momentum ground state is n = 2. For the non- interacting 
case of /i = 7r/2, the gap closes and n = N which correctly implies free acceleration in the 
thermodynamic limit. 

For fi > tt/2, the situation is more complicated: As has already been noted in there 
are special values of the interaction strength fJ^q = {Q — VjTi/Q at which the winding number 
becomes again macroscopic. In Fig. (|l|) we show as an example the spectrum at /is = 27r/3 
and = 12, M = 6. Note that not only the ground state gap closes at $ = 27r, but 
also various higher lying states become degenerate. Moreover, the closing of the gap at the 
special values /iq for Q = 2, 3, . . . , A^ occurs for all lattices sizes A^. 

In order to further understand this interesting behavior of the energy spectrum, we 
studied the Bethe Ansatz equations of the H-I chain with flux. The equations are given as 

i 

M 

N9{a,fi) = 27il + <!) + J20ia-a',fx), (2) 

a' 

where we have used Yang's change of variables for the pseudomomenta such that p = 
6{a,fi) = 2 arctan[cot(/i/2) tanh(a/2)]. The quantum numbers {/} are integers or half-odd 
integers and specify the states. The ground state is given by {/} = { — (M — l)/2, . . . , (M — 
l)/2}. 

Let us briefly recall the behavior of the roots {a} for the ground state as we adiabatically 
turn on the flux: The largest root Om goes to oo exactly at $ = 2(7r — /i) and the remaining 
M — 1 roots are distributed symmetrically around 0. The energy of the system is now equal 
to M — 1 particles on a chain of length A^ with $ = — the ground state in the = 1 sector. 
For < yU < 7r/2, further increase of $ renders om complex and it moves backwards along 
the line in with decreasing real part until at $ = 27r it has the value aM = i^c = ^i- Again, 
the remaining M — 1 real roots have been distributed symmetrically around 0. However, 
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the energy now is not related to the 5*^ = 1 ground state in a simple way. 

For 7r/2 < fi < 27r/3 and at $ = 2tt, we have two complex roots with zero real part, 
i.e. Ai^2 = ^(tt ± (tt — yu)), i.e. a two-string. And in general we find that at $ = 27r in the 
interaction interval {Q — 1)11/ Q < fi < {Q)'^/{Q + 1), we have a Q-string sitting on the 
imaginary axis symmetrically around m with individual roots given by = ivr + if3a = 
in + i{TT — fi){Q + 1 — 2a) + iSi- Sl represents an exponentially small finite-size correction, 
which is exactly zero for a 1-string, a 2-string and the two complex roots closest to iir in 
a Q-string with Q even. The remaining M — Q real roots are distributed symmetrically 
around 0. 

Note that at values of the interaction strength equal to fig, both a {Q — l)-string and 
Q-string coexist in the degenerate ground state. The energies associated with these strings 
are Eq^iiq) = and Eq^i^iiq) = 2(cos/iQ + 1) 0. This suggests, that for {Q - Vj-n/Q < 
fi < (Q)vr/ (Q + 1), when the Q-string corresponds to the ground state, the first exited state 
corresponds to the {Q + l)-string. 

We have therefore rewritten the BA equations (0) in order to incorporate the Q-strings 
sitting on the imaginary axis. The equations now are 

M-Q Q 

N9oo{a,fj.) =27il + J2 Ooo{a-a',2fx) + J20oiia,P',fJ^), (3a) 

a' 13' 
M-Q Q 

Neni/3,7T-fi)=0+ J2 ^io(«',A/i) + E^ii(/?-/^'>2(7r-/i)), (3b) 

a' 0' 

with 6'oo(a,/i) = 9{a,fi) the 

phase shift of real roots a, 9u{P,fi) = Re{2arctanh[tan(/3/2)/ tan(yu/2)]} the phase shift 
of complex roots and 6oi{a,P,fi) = Re{2 arctan[coth[(Q; — i[3)/2)/ tan(yu)]}, 6*10(0;, /?, /i) = 
— \ log[(cos(/? -|- 2/i) -|- cosh(a))/ (cos(/? — 2/i) -|- cosh(a))] the mixed phase shifts of a real with 
a complex root. The quantum numbers for the {/?}'s are zero and the fiux $ = 27r has been 
incorporated into the {/} quantum numbers. We caution the reader that care has to be 
taken in regard to the branch cuts of these phase shifts. The existence of roots on different 
branches is equal to a net shift of the quantum numbers and thus a different state. 
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In Fig. (0), we show the results of iterating the BA equations in the momentum 
71 sector for the ground state and the first excited state on a chain with M = 4 particles 
on = 8 sites. For < /i < 7r/2, the ground state corresponds to a 1-string state. As 
anticipated above, the first excited state is given by a 2-string state. At fi2, this 2-string state 
takes over as ground state and the 1-string state has vanished and in fact now corresponds 
to a higher lying state. An initially degenerate 3-string state appears and becomes the first 
excited state as /i is further increased. At /is, the 2-string now vanishes, the 3-string state 
takes over as ground state and a new 4-string state emerges. As there is no 5-string on a 
chain with M = 4 particles, this exchange of states ends here. 

We have plotted the energy difference AE = E'l — E'q at $ = 27r as a function of fi in 
Fig. (0). The closing of the gap at the special points fig can be clearly seen. Note that 
AE/Eq is already less than 2% at ;U = 77r/12 for this small system. For = 12 and M = 6, 
AE ~ 10"^ for the largest gap at /i ~ 77r/12 and decreases approximately exponentially as 

jj —>■ TV. 

As A^ oo, the gaps for ji ^ jiq should close as AE ~ A^""^. In Fig. (^, we have 
plotted 7 versus /i as extrapolated from calculations of up to A^ = 14 sites. Note that we 
can observe only Q-strings up to Q = 7 for these small sizes. Thus we can define 7 only 
for /i < = 67r/7. Moreover, for /i > /ig, useful data can come only from systems with 
M = Q + 1 particles. The result shows that 7 varies continuously for all < /i < vr. 
Alcaraz et al. have argued that both 1-string and 2-string for < /i < 7r/2 have the 
same scaling dimension. Therefore, the energy gap AE measures finite size behavior beyond 
simple conformal 1/A^ formulas and we interpret the continuous variation of 7 as indicating 
the presence of logarithmic corrections. A direct calculation of these corrections in the H-I 
chain has so far been done only at /i = and we are presently trying to extent these methods. 

However, at yU ~ 7r/2 (A ~ 0), we can calculate 7 directly by perturbation theory in A. 
As a starting point, we choose the plane wave basis of free particles at $ = and write 



We can then write the hopping term as 

N 27J- <|) 

T = -2 ^ cos(— m + —)dl^dm. (5) 

m=l 

Similarly, the interaction term is given as 

2A ^ /27r(m3 — m4)\ ^ ^ 

^ — 7^ cos I — I Ami+m3,m2+m4Wrrti'^m2'^»Ti3Um,4 5 

■^^ mi, 1712,1713,1714 \ / 

where the periodic Kronecker symbol is defined as An,m = 1 if n = m mod and otherwise. 

We now wish to identify the ground state \E'o aiid the first excited state \E'i which are 
degenerate at $ = 27r, i.e., 

^J)r^o|$=2. = ^It^i|$=2.. (7) 

Let M be odd and consider $ = 0. We then construct the ground state simply by filling 
all available momenta symmetrically, i.e., '^o = n!^~^(Af_i)/2 "^mlO)- Turning on the fiux, we 
find \E'o^\E'o|*=27r = \I'oT\E'o|*=o + 4sin(7r/A^). Keeping in mind Eqn. (0), we see that the first 
excited state is given by ^i|$=o = d^Sj^dcdD"^ o\<s>=o with A = -(M+3)/2, B = -(Af+l)/2, 
C = (M- l)/2 and D = (M-3)/2. A simple calculation shows that Eqn. holds indeed. 

For A 7^ 0, ^o|'i'=27r and \['i|$=27r will not be degenerate anymore, and the interaction V 
gives rise to a mixing term \E'gV\I'i. An explicit calculation shows that \E'oV^\E'i = ^[1 — 
cos27r/A^] ~ SAtc'^/N^. Therefore, we see that 7 = 3 at A = and $ = 27r. This agrees 
quite well with the numerical results of Fig. (|]) which were obtained for small systems of 
up to A^ = 14 sites and M = 7 particles |^. 

III. CLOSING OF ENERGY GAPS AT 1/r-FILLING 

Let us now study the behavior of the energy gap in the momentum vr sector away 
from half-filling. As mentioned above, one generally expects the winding number of the 
ground state to be = r for a 1/r filled band. As we have shown, this is true for the 
half-filled band except at isolated values of the interaction strength. For fillings equal to 
1/3, 1/4, 1/5, 1/6, .. ., we have explicitly calculated the energy spectrum as a function of fiux 



and we find that the gap closes at $ = vrr. In Fig. (H), we show a plot of the 1/3 filled 
case with = 77r/12 obtained by exact diagonalization of the H-I chain. This interaction 
strength corresponds to a large gap for the half- filled case. For 1/3 filling, however, the gap 
at $ = 37r is closed. Note that again various other degeneracies in the energy spectrum occur 
at this flux value. This behavior is unique to the 1/r filled chain — and by particle-hole 
symmetry to the (r — l)/r filled chain. It does not occur at all possible rational fillings such 
as, e.g., 2/5. 

The numerical results are again confirmed by perturbation theory. Let us, e.g., look 
at the 1/3 filled case with M odd and = 3M. The ground state \I^o is given as before 
and we need to identify the state which is degenerate with \l/o at $ = 37r. We find 
^il^^o = d\d\d}cdDdEdF^QW=() with A = -(M + 5)/2, B = -(M + 3)/2, C = -(M + l)/2, 
D = {M -l)/2, E = {M - 3)/2 and D = {M - 5)/2. 

For A 7^ 0, the exact diagonalization implies that \l/o|*=37r and \&i|$=37r will remain 
degenerate. Due to the structure of \I'i, we see that the mixing term indeed vanishes in first 
order in V. The second order mixing term is given by 

^¥,v ■ v^i = 

2A2 

-j^ cos27r(i3 - U)/N Ai^+i^^i^+i^ ■ cos27r(j3 - j4)/N Aj^+j^j^+j^ x 

il ,j2 ,j3 ,j4 

{'^oldl-^dj.^d'j^d'j.^di^di^dj^dj^dl^d^^dljdDdEdFl'^o)- (8) 

Explicitly performing the 576 possible contractions, we find that the diagrams arrange them- 
selves into 36 groups corresponding to a given ordering of contractions over the indices 
{A, B,C, D, E, F}. Only one index is left to be summed over and it is exactly this sum 
that gives zero for all the 36 groups individually. Therefore, we have an exact cancellation 
of second order diagrams for all fi in the critical region too. A generalization of this 
calculation to other 1/r fillings is straightforward, but becomes increasingly tedious as the 
number of possible contractions increases as {4 + r)\. 
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IV. CONCLUSIONS 



We have presented numerical and perturbative evidence for degeneracies in the momen- 
tum TT sector of the 1/r filled critical H-I model on a finite ring. Adding a non-integrable, 
say, next-nearest neighbor interaction W = X^^i to the Hamiltonian (|ID, we find 

that all the degeneracies discussed above for fi ^ 7i/2 are lifted and therefore all gap-closings 
disappear. The reader may now speculate that it is the integrability of the H-I chain which 
leads to the closed gaps. However, preliminary results from a related study of the repulsive 
Hubbard chain show that there the gap closes at f/ = only. It is an interesting ques- 
tion, if our results can be found also in other one-dimnesional quantum systems such as the 
long-ranged Haldane-Shastry chains. 

Furthermore, we have shown by exact diagonalization and iteration of the BA equations 
that the gap at $ = 27r (P = vr) and half-filling scales as power of the system size with 
variable exponent 7. An analytical calculation of 7 is in preperation. 

We close this paper by noting that as shown in M, the here considered structure of low- 



lying states in the H-I chain is qualitatively the same in the SC model |T0|. Therefore the 
results presented here will also hold in the SC model, up to the renormalization of quantities 
such as the spin wave velocity Vg. 
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FIGURES 

FIG. 1. The full spectrum of the H-I chain for A/" = 12 and M = 6 at /its = 27r/3. Note 

the level crossing for ground state and first excited state at <1> = 2tt. Various other level crossings 
enhance the periodicity of the ground state such that its winding number is n = 6. 

FIG. 2. Energy of the ground state and first excited state and their difference at $ = 27r 
(P = tt) for AT = 8 and M = 4. Note the closing of the gaps at 112 = 7r/2, /X3 = 27r/3 and 
/Lt4 = 37r/4. The finite gaps in the regions < A* < A^s and 1^3 < < i^a are just visible. 

FIG. 3. Assuming AE ^ N~'', we extrapolate 7 from finite size data of up to A?^ = 14 sites 
and half-filling. The error bars get larger as /i — ^uy since we need M > Q for jj, > jiQ-i- Data 
points for ii> ji'j = 67r/7 are therefore not expected to obey the simple power-law behavior. 

FIG. 4. The full spectrum for the H-I chain for AT = 12 and M = 4 at |Lt = 77r/12. Note 
the level crossing for ground state and first excited state at # = 37r (P = vr). Various other level 
crossings enhance the periodicity of the ground state such that its winding number is n = 6. 
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